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On Taub-NUT spacetime, we investigate hydrodynamic properties of perfect fluid spiraling inward 
toward the spacetime along a conical surface. On the equatorial plane of the Taub-NUT spacetime, 
we derive radial equations of motion with effective potentials and the Euler equation for steady state 
axisymmetric fluid. Higher dimensional global embeddings are constructed inside and outside the 
event horizons of the Taub-NUT spacetime. We also study the effective potentials of particles on 
the Taub-NUT spacetime in terms of gravitational magnetic monopole strength of the source, total 
energy and angular momentum per unit rest mass of the particle. 
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I. INTRODUCTION 



CO ' For many decades, four dimensional solutions of the Einstein equations have been extensively investigated in 

Cn I gravity community. Recently solutions of the Einstein equations or the coupled Einstein-matter equations in higher 

^_^ '■ dimensions have attracted many physical interests since the discovery of supergravity and superstring theory. Taub- 

Q ■ NUT metric flJ] is a local analytic solution of the vacuum Einstein equations. When the metric is expressed in 

O^' Schwarzschild-like coordinates, we have coordinate singularity occurring at certain values of the radial coordinate 

; I ' where g^r becomes infinite and corresponds to bifurcate Killing horizons. The Taub-NUT spacetime is involved in 

bJO, many modern studies in general relativity. Hawking has proposed that the Euclidean Taub-NUT metric could give rise 

'~~' ' to the gravitational analogue of the Yang-Mills instanton [3] . In this case the Einstein equations are fulfilled with zero 

,__! ' cosmological constant and the manifold is R^ with a boundary which is a twisted three-sphere possessing a distorted 

^ ' metric. The Kaluza-Klein monopole has been obtained by embedding the Taub-NUT gravitational instanton into 

^^ five dimensional Kaluza-Klein theory [3, |j| . The non-Abelian target space duals of the Taub-NUT spacetime have 

^^ I been also studied in terms of the local isometry group SU(2)xU(l) [3]. The Taub-NUT spacetime has been shown 

^~~^ ' to be related with SU(2) calons through T-duality [S]. Carter has shown that the Hamilton- Jacobi equation for the 

'^ , geodesies in the Taub-NUT metric separates in certain coordinate systems [3] . The gravitomagnetic monopole source 

lO ' effects have been also studied in the Taub-NUT spacetime Q. In the Kerr-Taub-NUT-de Sitter metrics, separability of 

^^ ' the Hamiltonian- Jacobi equation has been studied in higher dimensions 9] . Recently, a rotating Schwarzschild black 

hole has been studied to investigate effective potentials for null and timelike geodesies of particles and hydrodynamics 

associated with general relativistic Euler equation for the steady state axisymmetric fluid [lOj . 

A familiar feature of exact solutions to the field equations of general relativity is the presence of singularities. 
As novel ways of removing the coordinate singularities, the higher dimensional global flat embeddings of the black 
hole solutions are subjects of great interest both to mathematicians and to physicists. It has been well-known in 
J^ ' differential geometry that four dimensional Schwarzschild metric [11] is not embedded in R^ [ij]. Recently, (5+1) 
dimensional global embedding Minkowski space (GEMS) structure for the Schwarzschild black hole has been ob- 
tained [ij] to investigate a thermal Hawking effect on a curved manifold [IJl associated with an Unruh effect [l^ 
in these higher dimensional spacetime. In (3-Kl) dimensions, the global flat embeddings inside and outside of event 
horizons of Schwarzschild and Reissner-Nordstrom black holes, have been constructed and on these overall patches 
of the curved manifolds four accelerations and Hawking temperatures have been evaluated by introducing relevant 
Killing vector s 11611 . Recently, the GEMS scheme has been applied to stationary motions in spherically symmet- 
ric spacetime |l3], and the Banados-Teitelboim-Zanelli black hole [T^ has been embedded in (34-2) dimensions to 
investigate the SO(3,2) global and Sp(2) local symmetries [l9| . 

In this paper, exploiting the fact that the coordinates t and are cyclic in the Taub-NUT spacetime we find 
the timelike Killing field and the axial Killing field, to which we can obtain the conserved energy and the angular 
momentum per unit rest mass for geodesies. In the Taub-NUT spacetime we then investigate hydrodynamic properties 
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of the perfect fluid spiraling inward toward the Taub-NUT spacetime along a conical surface. On the equatorial plane 
of the Taub-NUT spacetime, we derive the radial equations of motion with the effective potential. We also study the 
effective potentials of particles on the Taub-NUT manifold in terms of gravitational magnetic monopole strength of 
the source, total energy and angular momentum per unit rest mass of the particle. Finally, we construct the higher 
dimensional global embeddings inside and outside the event horizons of the Taub-NUT manifold. 

This paper is organized as follows. In section II we introduce the Taub-NUT spacetime to study its hydrodynamics, 
and in section III we construct the global flat embeddings inside and outside of event horizon of the Taub-NUT 
spacetime. Section IV includes summary and discussions. 

II. HYDRODYNAMICS OF TAUB-NUT SPACETIME 

We start with the Taub-NUT 4-metric given by 

A E 

ds^ = (dt + 21 cos edcj))^ + —dr"^ + Y-ide"^ + sin^ ed(j)'^) (2.1) 

where, in the units G = c = 1, 

A = r^ -2Mr-l'^, T,^r'^ + P, (2.2) 

with the parameters M and I being associated with the mass and gravitational magnetic monopole strength of the 
source. Defining r± as 

r±^M±{M^ + ff/^ (2.3) 

we can rewrite A as 

A = (r-r+)(r-r_). (2.4) 

Here note that r+ > / is an event horizon but r_ is not since r_ is negative. 
The four velocity is given by 

u^ = ^> (2-5) 

where we can choose r to be the proper time (affine parameter) for timelike (null) geodesies. From the equation of 
motion of a test particle in the Taub-NUT spacetime, the particle initially at rest at infinity spiral inward toward 
the spacetime along a conical surface of constant 9 = 6^0 where 6^0 is the polar angle at infinity. For a fluid which 
is at rest at infinity and approaches supersonically to the spacetime, one may take the approximation to simplify the 
hydrodynamical equations 

u^ = ^ « 0. (2.6) 

dr 

As in the Schwarzschild black hole since the coordinates t and are cyclic we have the timelike Killing field ^"^ and 
the axial Killing field V^". Corresponding to the Killing fields ^° and -0" we can then find the conserved energy E and 
the angular momentum L per unit rest mass for geodesies given as follows 

E = -gabCy" 

{r-r+){r-r^) t , „, 0. 
= ^a^2 (^ +2/7.*), 



2l{r-r+){r-r^) ^, , /, 3 , ,2^ .,2 n 4?2(r - r+)(r - r_) 



,2 + p "* + W' + n Bin^ e - ^ ^^Y;^ ^ cos^ ] u\ (2.7) 



where m° are four velocity of the locally nonrotating observers defined by (j2.5p . Moreover, we can introduce a new 
conserved parameter k defined as 



abu'^u'' (2.8) 



whose values are given by k = 1 for timelike geodesies and k = for null geodesies. 

In the case of geodesies on the equatorial plane 9 — 7r/2, m* and u'^ are given in terms of E and L as follows 



»' 



( (r-'^Cf ^_^,^ 



(r^ + P)2 + AP{r — r^){r — r^) \ (r — ?'+)(7' — f-) 

^^"^ = r^5 Wr^ 7777 77 -AL + 21E), (2.9) 

which are substituted into (j2.8p to yield the radial equation for the particle on the equatorial plane 

]^E^ ^^^u'-u^ + V{r,E,L), (2.10) 

with the effective potential 

2(r2 + P) ^2[(r2 + Z2)2+4/2(r~r+)(r-r_)]^ ^ ^ ■ ^^^ 

For the null geodesies with k = 0, the effective potential V[x,y) for the particles with the total energy per unit 
rest mass E = Q.l and angular momentum per unit rest mass L = 0.2M is shown in the first graph in Fig. 1 where 
X and y denote the dimensionless variables x = r/M and y = l/M, respectively. We observe here that there exist 
many peaks on the potential surface. In the third graph, we plot the surface z = (x^ + y^)^ + 4y^(a;^ — 2x — y"^) and 
the plane z = 0, which shows the intersection curve corresponding to the vanishing values of the denominator of the 
second term in (J2.11I) . Along this curve we have peaks of the potential surface in the first graph. In the second graph 
in Fig. 1, V{x,y) is shown for the timelike geodesies of the particles with E = 0.1 and L — 0.2A/. In this case there 
exist peaks similar to those in the null geodesies. The potential surface is deformed due to the first term proportional 
to K, in (J2.11I) , which however does not affect the peak curve in the third graph since the denominator of this term is 
positive definite. 

The fundamental equations of relativistic fluid dynamics can be obtained from the conservation of particle number 
and energy-momentum fluxes. In order to derive an equation for the conservation of particle numbers one can use 
the particle flux four vector nu° [2^ 

Va(nu'^) = ^=Va(\/^ nu") ^ 0, (2.12) 



where n is the proper number density of particles measured in the rest frame of the fluid and Vq is the covariant 
derivative in the Taub-NUT curved manifold and g = det gab- For steady state axisymmetric flow, the conservation 
of energy-momentum fluxes is similarly described by the Einstein equation [2l| 

^bTl = ^Vb(V^ Tl) = 0, (2.13) 
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where the stress-energy tensor T°'^ for perfect fluid is given by 

T""^ = pu'^u'' + ig"'' + u''u'')P (2.14) 

with p and P being the proper internal energy density, including the rest mass energy, and the isotropic gas pressure, 
respectively. The Einstein equation ()2.13p can be rewritten in another covariant form 

Ua^biiP + P)u'') + {P + p)u^VbUa + V^P = 0. (2.15) 

On the other hand, from (|2.ip . (|2.2p and (|2.3|) we can find the inverse of the Taub-NUT metric 

r2+/2 AP{r^+P) „ {r-r^){r-rJ) ,, 1 



{r~r+)(r-rA [(r^ -^ F)2 + 4?2(r - r+)(r - r_)] sin^ 0' ^ r^ + P ' '' r^+P' 

,.. . !^±l ,*0 . ^^(-^ + ^^) (2 16) 

y [{r^+py+iP(^r-r+){r-rA]sm^e' ^ [{r^ + P)^ + 4P{r - r+){r - rA]sm^ 9' ^' ' 

Exploiting the inverse 4-metric g"'^ in (J2.16I) and acting the projection operator gab + UaUb on the equation (I2.15P we 
can obtain the general relativistic Euler equation on the direction perpendicular to the four velocity 

(P + p)u^VbUa + {gab + UaUb)V^P = 0. (2.17) 



After some algebra, from (|2.17p we obtain the radial component of the Euler equation for the steady state axisymmetric 
fluid 

d ^ ^ 2r — r+ — r_ 2(r — r+)(r — r_)r ^ APB 



dr' ' r^+P (r2 + Z2)2 (r2+/2)^ 



(r — r+)(r — r^)A P + p \ r^ + P J dr 

where 

A ^ (r^ +/2)2sin2 6*- 4/2(r-r+)(r-r_) cos^ 61 

^ (r2 + P)(r - r+)(r - r_)[(2r - r+ - r_)(r2 + P) -- 4r(r - r+)(r - r_)] 

(r2 + Z2)2+4p(^_j,^)(^_^_) ■ 1^ • ^^ 

In the vanishing I limit, the result (12.18^ is reduced to that of the (rotaing) Schwarzschild black hole |l(j |. 
Multiplying p.lSp by m° we can project it on the direction of the four velocity to obtain 

nu- (Va (^^) - ll^aP] = 0, (2.20) 

where the continuity equation (|2.12p has been used. The radial component of (|2.20p yields 

dp P + pdn A - r 



dr n dr 



(2.21) 



Here the energy loss A and the energy gain T are introduced to set the decrease in the entropy of inflowing gas equal 
to the difference A — F. Here it is interesting to see that the result (j2.2ip holds also in the (rotating) Schwarzschild 
black hole tlO|]. 

III. GLOBAL EMBEDDINGS OF TAUB-NUT SPACETIME 

After tedious algebra, for the Taub-NUT spacetime in the region r > r+ we can obtain the (6+5) global embedding 
Minkowski space (GEMS) structure 

ds^ = -{dzy + [dz^f - {dz^f + [dz^f + {dz^f + {dz^f + {dz^'f - [dz'f - {dz'f + [dzy - {dz'°f (3.1) 

with the coordinate transformations 

^0 ^ ^_i /^r-r+)(r-r_2\ ^^s - sinh«:(t + 2/0), 
\ r^ + i^ / 2 

-^ / (r — r_|_)(r — r_) ^ '^ 

r2 + /2 J -™2 

1/2 



'^"^ ( ^ T]2 I ^'^^ o cosh K{t + 2l(j)), 



1/2 



z = K. ( ^ -^ j sin -sinhK(t — 2Z0), 



z — K T, -t; sin - cosh K(i — 2/(/)), 

* r2 + /2 / 2 ^ ' 



— K 

4 ,'(r2+/2)2+4/2/^_ W^_^_)xl/2 

z^ = I ^ „ ,„ 1 sin«cos( 

r^ + P 

5 / (r2 + /2)2 + 4;2(^ _ )(^ _ ^_) X 1/2 

z^ = I -^^ 77-^ — ^ '^ 1 smfsmc 

r2 +P 



z 



6 



(r2 + /2)2 _^ 4^2(^ - r+)(r - r-) ^ ^^^ 

r2 +;2 



cos 6*, 



,^^i4,^+VA"7"--t"';:"-')'^o.<', 



4 / V r2 + P 



z 



8 



4/2 + t^-^ 
4 



1/2 



(r ~ r j^)(r — r J) 

y,2 _|_;2 



1/2 



sin0, 



-,9 



{p(r)[(r+ - r^){r'^ + /2)2 + (r^ + l-^)s{r)\ + g(r)(r^ + /2)(r+ + r,)/^}!/^ _ 

(r - r_)i/2(r_^ _ r_)(r2 + Z2)3/2 " AH, 

10 _ / ^ /p(r)(r^+/2)(r++r_)/2 + 5(r)[(r+-r_)(r2+/2)2 + (^2_^;2)g(^)] 



z" — dr 
z^" = dr 



(r — r^){r^ — r^)^(r^ + P)^ 
[r{P + lY + 2/2(2r - r+ - r_)(r2 + f) - APrjr - r+)ir - r,) ]^ ^^ ^^^ 

[(r2 + /2)2 + 4^2(^ _ ^_^)(^ _ r_)](r2 + ;2)3 



where k is given by 



9ir), (3.2) 



16 



-^+^ ^ +,2 



(3.3) 



and the positive definite functions are given by 



p{r) = (r+ — 7'_)r +?'+(?'+ — r_)r +/ r+r-, 
q{r) = (2r^ + 3/2)r_r + /2(r;j_+r+r_ +2/2), 

s(r) = {r-^- — r^)r +2{l — r-^-r^)r. 



(3.4) 



Here note that k in p.3p is not the surface gravity, different from those in the embeddings in the static black holes 
without the shift functions [H, [H, [H, [Hi • 

Next, in order to investigate the region r < r+ we rewrite the Taub-NUT 4-metric as 



ds^ = ( ] {dt + 21 cos ed(t)f- (-^) dr"^ + T.{de^ + sin^ edcj)^) 

in terms of the positive definite lapse function inside the event horizon r_|_ 

A (r+ — r){r — r-) 
" S " r2 + /2 

to yield the (6+5) GEMS structure p. II) with the coordinate transformation 

^0 ^ ^-1 ( (r+^ZliL^lA] cos^cosh«(t + 2/0), 
^1 ^ ^-1 f v'+-;a-j--; I ^Qg ^ gj^j^ ^(^ ^ 2/^), 

^2 ^ ^-if ^'+~:^^'.7''' ] sin^cosh«;(i-2/0), 



V 


r2 + /2 






-i(ir+ 


-r)(r- 


r- 


-) 


\ 


r2 + /2 






-ifir+ 


-r)(r- 


r- 


-) 


\ 


r2 + /2 






-ifir+ 


-r)(r- 


r- 


-) 
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1 



r2 + / 

1/2 



1/2 



sin — sinh K{t — 2/0), 



Z^ = 4/2 + -^-2 



(r+ — r)(r — r_) 



1/2 



^8 = 4/2 + ±^-2 



4 7 V r2 + /2 

4" j I r2 + /2 



COS! 



smty, 



z' = fir), 
zi" = 5(r), 



(3.5) 



(3.6) 



(3.7) 



where (z , z , z ), /(r) and g{r) are the same as those in 



IV. CONCLUSIONS 

In conclusion, we studied the Taub-NUT spacetime to investigate hydrodynamic properties of the perfect fluid 
whirling inward toward the Taub-NUT spacetime along a conical surface. On the Taub-NUT manifolds we then 
constructed the (6-1-5) dimensional global embeddings inside and outside the event horizons of the manifold. We 
also studied the effective potentials of particles on the equatorial planes of the Taub-NUT spacetime in terms of the 
gravitational magnetic monopole strength of the source, total energy and angular momentum per unit rest mass of 
the particle. 
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